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Abstract 

A polytope D whose vertices belong to a lattice of rank d is Delaunay if there is a 
circumscribing d-dimensional ellipsoid, £, with interior free of lattice points so that the 
vertices of D lie on £. If in addition, the ellipsoid £ is uniquely determined by D, we 
call D perfect. That is, a perfect Delaunay polytope is a lattice polytope with a circum- 
scribing empty ellipsoid £, where the quadratic surface <9£ both contains the vertices 
of D and is determined by them. We have been able to construct infinite sequences of 
perfect Delaunay polytopes, one perfect polytope in each successive dimension starting 
at some initial dimension; we have been able to construct an infinite number of such 
infinite sequences. Perfect Delaunay polytopes play an important role in the theory of 
Delaunay polytopes, and in Voronoi's theory of lattice types. 

1 Introduction 

Consider the lattice and a lattice polytope D. Then, D is a Delaunay polytope if it can be 
circumscribed by an ellipsoid 8 so that: (1) £ has no interior Z d -elements, and (2) the only 
Z d -elements on £ are the vertices of D. We will say that the circumscribing ellipsoid £, with 
these two properties, is an empty ellipsoid. If D is a Delaunay polytope with circumscribing 
empty ellipsoid £, then D is perfect if the quadratic surface 88 is uniquely determined by 
the vertex set of D. 

If D is Delaunay, with empty ellipsoid £, the vertices of D are given by 88 f] D = 88 f] Z d . 
Typically, the quadratic surface 88, through this vertex set, can be deformed continuously. 
It is only in the case where D is a perfect Delaunay polytope that this surface is uniquely 
determined by 88 f]D. We refer to 88 as a perfect ellipsoid, and by extension, to refer to 
D perfect Delaunay polytope. 

We have studied perfect Delaunay polytopes by constructing infinite sequences of them, 
one perfect Delaunay polytope in each successive dimension, starting at some initial dimen- 
sion. We have been able to construct an infinite number of infinite sequences of perfect 
Delaunay polytopes. 

The perfect Delaunay polytopes in each sequence have similar combinatorial and geo- 
metric properties. One of our constructions is the sequence of G 6 -topes, Gf, d = 6,7,... , 
with the initial term the semiregular 6-dimensional Gossett polytope G 6 with 27 vertices 
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(2 2 i in Coxeter's notation); each GVtope is asymmetric with respect to central inversion, 
and Gq has ^ttli 4- d vertices. Each term in the sequence is a combinatorial analogue of 
G% = 221. Another of our constructions is the sequence of Gf-topes, G%, d = 7, 8, with 
the initial term the 7- dimensional Gossett polytope G7 with 56 vertices (331 in Coxeter's no- 
tation); each GV-tope is centrally-symmetric, and Gj has ^( d 1 2 1 ) vertices. Each term is again 
a combinatorial analogue of the initial Gosset polytope Gj = 3 3 i. Just as the 6-dimensional 
Gossett polytope can be represented as a section of the 7-dimensional one, each Gq of the 
asymmetric sequence can be represented as a section of the G d+1 of the symmetric sequence. 

Properties and Connections: Perfect Delaunay polytopes are fascinating geometric ob- 
jects. There is a number of other reasons why we studied these polytopes. 

1. Perfect Delaunay polytopes are maximal Delaunay polytopes - it follows immediately 
from the definition that a perfect Delaunay polytope cannot fit properly inside another 
Delaunay polytope of the same dimension. 

2. Each perfect Delaunay polytope D uniquely determines a positive definite quadratic 
form on Z d , which in turn, determines a Delaunay tiling of Z d ® R with D as one of its 
Delaunay tiles (details on how these tilings are constructed are given in the following 
section). Therefore, each of our sequences of perfect Delaunay polytopes determines 
a corresponding sequence of Delaunay tilings. Just as the perfect Delaunay polytopes 
in the sequence are combinatorial analogues, the corresponding Delaunay tilings have 
somewhat similar combinatorial properties. 

Each perfect Delaunay polytope in the sequence of G 6 -topes uniquely determines a 
Delaunay tiling that is an analogue of that for the root lattice E % in E 6 . Similarly, 
each term in the sequence of G 7 -topes uniquely determines a Delaunay tiling that is 
an analogue of that for the root lattice E 7 in E 7 . 

3. A quadratic form ip in d indeterminates is called perfect (after Voronoi) if it is positive 
definite and is uniquely determined by its minimal vectors through the linear system 
of equations {x[m] = min^oV 9 I m is a minimal vector for <p}, where the coefficients 
of x are the variables, min^o V 9 i s the minimal value achieved by (p on the non-zero 
elements of Z d , and each equation corresponds a minimal vector m G Z d for (p. If (p is 
perfect, then the only solution is x — ¥■ 

There is an inhomogeneous analogue to perfect form that plays an important role in 
the study of Delaunay polytopes. Let £ be an empty ellipsoid around a Delaunay 
polytope D, defined by inequality /(x) < 0. Since £ is empty, / is non-negative on Z d , 
and assumes its minimum value zero on <9£ f] Z d . If D is perfect then the coefficients 
of the inhomogeneous quadratic function / are determined, up to a scale factor, by 
the linear system of equations {/(z) = | z e 8,f]Z d }. The rank of this linear 
system must then be ^tbl) _|_ ^ which requires a perfect Delaunay polytope to have 
at least this number of vertices. The vertices of a perfect Delaunay polytope are used 
to determine the coefficients of / in much the same way that the minimal vectors are 
used to determine the coefficients of a perfect form. Therefore, when D is perfect it 
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is natural to refer to / as a perfect quadratic function, and to refer to £ as a perfect 
ellipsoid. Clearly is perfect Delaunay in Z° and [0, 1] is perfect Delaunay in Z 1 . Next 
perfect Delaunay polytope lives in dimension 6. 

The conditions that a perfect quadratic function must satisfy are more severe than 
those for a perfect form, which requires at least n (" 2 +1 ) pairs of minimal vectors. This 
leads to the speculation that the growth of numbers of arithmetic types of perfect 
quadratic functions with dimension will be much slower than for numbers of types of 
perfect forms. This has been borne out by the data now available on perfect Delaunay 
polytopes (see http : //www, liga. ens . f r/~dutour/p . In fact, the growth of numbers 
of arithmetic types of perfect function seems much slower that that for most classes of 
geometric objects studied in geometry of numbers. 

Perfect forms have been introduced by Korkine and Zolotareff (1873), who classified all 
the perfect forms through dimension five (the term perfect was coined by Voronoi). In 
addition, they found several perfect forms in higher dimensions, including the forms cor- 
responding to the root lattices A n , D n (both for n EN), E 6 , E 7 , and E s , although they 
used completely different roman letters to denote these lattices. Conway and Sloane 
(1988) give a survey of perfect forms through dimension seven. More recently, the study 
of perfect forms has been taken up by Martinet who published a comprehensive mono- 
graph on the topic in 2003. As of the end of 2006 perfect forms have been classified in 



dimensions through eight (see http://fma2.math.uni-magdeburg.de/~latgeo/) 



4. A significant role for perfect Delaunay polytopes is to provide geometric labels for a 
class of edge forms for the homogeneous domains of Delaunay tilings of 7L d ® R. These 
domains, also referred to as L-type domains, are polyhedral cones of positive forms, 
and the edge forms are those that lie on the extreme rays. L-type domains have been 
first introduced and studied by Voronoi, and they play a central role in his classification 
theory for lattices, his theory of lattice types. Delaunay associated these polyhedral 
cones with Delaunay tilings. Prior to the discovery of the sequence {Gf} of GVtopes 
by Erdahl and Rybnikov in 2001 only finitely many arithmetic types of edge forms 
of full rank were known. All of the infinite sequences of perfect Delaunay polytopes 
that we report on correspond to infinite sequences of edge forms of full rank for L- 
type domains. Later infinite series of perfect Delaunay polytopes were constructed 
by Dutour (2005) and Grishukhin (2006). Our multiparametric series overlap with 
the 1-parametric series of Dutour only in small dimensions (6,7,8). Our series contain 
all of the series of Grishukhin, although his approach is different. The remarkable 
property of Dutour's series is that the size of the vertex set of Dutour's polytope grows 
exponentially with the dimension, for his (i-polytope ED^, for even ef s, contains a 
section isometric to the (d — l)-halfcube (for odd d's it has 2 sections isometric to the 
(d — 2)-halfcube). Dutour conjectured that this polytope ED^ has the largest number 
of vertices among all perfect Delaunay polytopes of dimension d. 

Perfect Delaunay polytopes, and the infinite cylinders with perfect Delaunay polytopes 
as bases, also provide labels for edge functions for the inhomogeneous domains of 
Delaunay polytopes. Associated with each lattice Delaunay polytope is a polyhedral 
cone of inhomogeneous quadratic functions, and we refer to the functions lying on 
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extreme rays as edge functions. These domains play an important role in the structure 
theory of Delaunay polytopes. 

5. Delaunay tilings of lattices and Voronoi's theory of L-types are intrinsically connected 
to the theory of moduli spaces of abelian varieties. On a philosophical level this 
connection stems from the idea of parametric approach to the study of geometric 
objects. In geometry of numbers and number theory this approach was pioneered by 
Hermite (1850). Korkine, Zolotareff, and Voronoi referred to this approach as Hermite's 
method of continuous variation of parameters. It is worth noting that Voronoi's interest 
in geometry of numbers developed as a result of his work on irrationalities of the third 
degrees and, in particular, elliptic curves. The curous reader may consult Delone 
[Delaunay] and Faddeev (1964). For recently discovered connections between moduli 
spaces of abelian varieties and Voronoi's theories of perfect and of L-types see Alexeev 
(2002) and Shepard-Barron (2005). We are not aware for any results interpreting 
perfect perfect Delaunay polytopes and L-types defined by them, but we certainly 
anticipate them. 

The following section is devoted to the details of how perfect Delaunay polytopes 
relate to these fundamental polyhedral cones in geometry of numbers, the homogeneous 
domain of a Delaunay tiling, and the inhomogeneous domain of a Delaunay polytope. 

Historical notes: Perfect Delaunay polytopes were first considered by Erdahl (1975) in 
connection with lattice polytopes arising from the quantum mechanics of many electrons. 
He observed (1975) that Delaunay tilings of 0- and 1-dimensional lattices consist entirely 
of perfect Delaunay polytopes, and showed that the Gosset polytope G 6 = 2 21 with 27 
vertices was a perfect Delaunay polytope in the root lattice E$. He also showed that there 
were no perfect Delaunay polytopes in dimensions 2, 3, and 4. These results were further 
extended by Erdahl (1992) by showing that the 7-dimensional Gosset polytope Gj = 3%i 
with 56 vertices is perfect, and that there are no perfect Delaunay polytopes of dimension d 
for 1 < d < 6. Erdahl also proved that Gq and G7 are the only perfect Delaunay polytopes 
in the Delaunay tilings for the root lattices. Deza, Grishukhin, and Laurent (1992-1997, in 
various combinations) found more examples of perfect Delaunay polytopes in dimensions 15, 
16, 22, and 23, but all of those seemed to be sporadic. 

The first construction of infinite sequences of perfect Delaunay polytopes was described 
at the Conference dedicated to the Seventieth Birthday of Sergei Ryshkov (Erdahl, 2001), 
and later reported by Rybnikov (2001) and Erdahl and Rybnikov (2002). Perfect De- 
launay polytopes have been classified up to dimension 7 - Dutour (2004) proved that 
Gq = 2 2 i is the only perfect polytope for d — 6. It is strongly suspected that the ex- 
isting lists of seven and eight dimensional perfect Delaunay polytopes are complete (see 
http : //www . liga . ens . f r/~dutour ) . 



2 Quadratic Polynomials on Lattices 

Perfect Delaunay polytopes play an important role in the L-type reduction theory of Voronoi 
and Delaunay for positive (homogeneous) quadratic forms. An L-type domain is the col- 
lection of all possible positive quadratic forms that give the same Delaunay tiling D for 
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Z d . L-type domains are relatively open polyhedral cones, with boundary cells that are also 
L-type domains - these conical cells fit together to tile the cone of positive quadratic forms, 
which is described in the next section in detail. Simplicial Delaunay tilings label the full 
dimensional conic "tiles", and all other possible Delaunay tilings label the lower dimensional 
cones. 

The significance of extreme L-types is much due to their relation to the structure of 
Delaunay and Voronoi tilings for lattices. The Delaunay tilings that correspond to edge 
forms play an important role: The Delaunay tiling for an L-type domain is the intersection 
of the Delaunay tilings for edge forms for the L-type (Erdahl, 2000). There is a corresponding 
dual statement on the structure of Voronoi polytopes: The Voronoi polytope V v for a form <p 
contained in an L-type domain L is a weighted Minkowski sum of linear transforms of Voronoi 
polytopes for each of the edge forms of £ . The latter dual result was first established by 
H.-F. Loesch in his 1990 doctoral dissertation, although it was first published by Ryshkov 
(1998, 1999), who independently rediscovered Loesch's theorem; this dual result was given 
a shorter and simpler proof by Erdahl (2000). 

Edge forms that are interior to the cone of positive forms are rare in low dimensions. 
They first occur in dimension 4: D 4 has an extreme L-type. A good proportion, but not 
all, of the edge forms appearing in lower dimensions relate either directly or indirectly to 
perfect Delaunay polytopes. As shown by Dutour and Vallentin (2005) this situation does 
not persist in higher dimensions: there is an "explosion" of the number of inequivalent 
interior edge forms in six dimensions, and only a tiny fraction of these are inherited from 
perfect Delaunay polytopes. 

The Voronoi and Delaunay tilings for point lattices are constructed using the Euclidean 
metric, but are most effectively studied by injectively mapping the lattice into Z d , and 
replacing the Euclidean metric by an equivalent metrical form. For a d- dimensional point 
lattice A with a basis b 1; b 2 , h d this is done as follows. A lattice vector with coordinates 
zi, Z2, Zd relative to this basis can be written as z = Bz, where B = [bi, b 2 , bd] is the 
basis matrix and z is the column vector given by [z\, Z2, zaY- The squared Euclidean length 
is given by |v| 2 = z*B*Bz =</?b(z). This squared length can equally well be interpreted as 
the squared length of the integer vector z G Z d relative to the form </? B - Therefore, the 
Dirichlet-Voronoi and Delaunay tilings for A, constructed using the Euclidean metric, can 
be studied using the corresponding Dirichlet-Voronoi and Delaunay tilings for Z d constructed 
using the form </?b as metric. Moreover, variation of the Dirichlet-Voronoi and Delaunay 
tilings for A in response to variation of the lattice basis can be studied by varying the metric 
ip for the fixed lattice Z d . In this section we will keep the lattice fixed at Z d , and vary 
the quadratic form <p. We call a form <p positive if <p[z\ > for any z, and we call it 
positive definite if tp[z] > for any z ^ 0. Same terminology is applied to all numbers and 
functions. Positive, but not positive definite forms are referred to as positive semidefinite, 
or just semidefinite. 

The inhomogeneous domain of a Delaunay polytope: Let T d be the cone of real 
quadratic polynomials defined by: 

9 d = {feR[x 1 ,...,x d ] | deg/ = 2, Vz G Z d : f(z) > }. 

The condition Vz G Z d : /(z) > requires the quadratic part of / to be positive, and requires 
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any subset of R d where / assumes negative values to be free of Z d -elements. The real quadric 
determined by the equation /(x) = 0, where / G 7 d , might be empty set or it might have 
the form 

88. f — <9£ x K, 

where £o is an ellipsoid of dimension n > and K a complementary subspace, of dimension 
d — n. The latter case is the interesting one - depending on the dimension of K, £j is either 
an empty (of lattice points) ellipsoid or an empty cylinder with ellipsoidal base. 

For any function / G 7 d we denote by V(/) the set { z G Z rf | /(z) = 0}. In the case 
where the surface <9£/ = {xGZ d Cg>R| /(x) = 0} includes integer points and is bounded, 
V(/) = £/ H Z d is the vertex set for the corresponding Delaunay polytope Df — convV(/). 
Conversely, if D is a Delaunay polytope in Z d ® R, there is a circumscribing empty ellipsoid 
determined by a function f D G 7 d . More precisely, since D is Delaunay, there is a form </?d, 
a center c, and a radius R, so that /d(x) = </?d[x — Cd] — < defines a circumscribing 
empty ellipsoid. Since /d is positive on Z d , it is an element of 7 d . 

In the case where £ / is an empty cylinder, there is an infinite number of integer points 
lying on this surface. When this happens V(/) is the set of all lattice points on the boundary 
of a non-bounded polyhedron Df — convV(/), which we call a Delaunay polyhedron (we 
follow a recent convention that a polyhedron deserves to be called a a polytope when it is 
bounded). Although elements ofV(f) need not be vertices in the sense of real polyhedral 
geometry, we still refer to them as such and denote them by vertDf. With this definition 
of Delaunay polyhedron, for a given form <p > Delaunay polyhedra need not form a tiling 
of Z d ® R. They do form a face-to-face tiling if and only if the real kernel of the quadratic 
part of / intersects with Z d over a sublattice of the same dimension, i.e. if dim &c? K (p = 
dim{^et]R <p n Z d }. Some authors refer to this conditions as that tp has a rational kernel or as 
that ip has a rational radical, although both expressions are somewhat ambiguous - a more 
accurate way of expressing this condition is to say that <p has fully rational kernel. 

Definition 1. A polyhedron D C Z rf ®Ris called a Delaunay polyhedron, if there is / G ^ 
such that D = convV(/). 



Definition 2. Let D be a Delaunay polyhedron in Z d ® R. Then, the (inhomogeneous) 

y d D = { f e T d | Df — D.} 



domain y d D for D is: 



Such domains are relatively open convex cones that partition y d . 

Each element / G < 3 >d ~ ) satisfies the linear equations /(z) = 0, z G D n Z d . When D is 
a single lattice point and d > 0, 3"^ is a relatively open cone of dimension ( d ~^ 2 )- When D 
is a 1-polytope, y d D is a relatively open facet of the partition with dimension dim T d — 1 = 
( d + 2 ) - 1. When the rank of the system {/(z) = | z G D n Z d } is full, i.e. ( d + 2 ), the only 
possible / is and D = Z d <g> R. If D is a rf-dimensional (d > 0) Delaunay simplex then 
dim^ = but if D is a perfect Delaunay polyhedron and d > 0, then dimJ"^ = 1. 

Definition 3. A function p G T d is perfect if the system of equations 



/(z) = minp 



z G V(p) 
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has a unique solution / = p. In this case we also call conv V(p) a perfect polyhedron. 

The elements of 0- and 1-dimensional inhomogeneous domains are perfect, and the De- 
launay polyhedra that determine such domains are perfect. The domain that consists of 
strictly positive constant functions corresponds to the empty Delaunay polytope. The per- 
fect subsets V(p) must be maximal among the subsets { V(/) | / G T d }. With the exception 
of the 0-dimensional perfect domains, all inhomogeneous domains J"^ have proper faces that 
are inhomogeneous domains of lesser dimensions. Perfect quadratic functions are analogues 
of perfect forms - both achieve their arithmetic minimum at a sufficient number of points so 
that the representations of the minimum uniquely determine the polynomial. 

The following theorem and the next subsection show the important role played by perfect 
Delaunay polytopes in the theory of lattice Delaunay polytopes and Voronoi's L-types. 

Theorem 4. (Erdahl, 1992) Let D be a Delaunay polyhedron in 7L d ® R. Then 

I {p | p is perfect, vert DCV(p)} 

This theorem shows that an arbitrary element / G J*d has the following representation: 

{p | p is perfect, vertDCV(p)} 

where Uf jP > 0. 

The homogeneous domain of a Delaunay tiling: Voronoi's classification theory for 
lattices, his theory of lattice types (L-types), was formulated using positive definite forms and 
the fixed lattice Z '. In this theory two lattices are considered to be the same type if their 
Delaunay tilings are affinely equivalent. This characterization is not apparent in Voronoi's 
memoirs and was brought to light by Delaunay. Here is a brief reminder of Delaunay's 
approach to Voronoi's theory of L-types. Consider a positive definite quadratic form ip. 
Then a Z d -polytope D is Delaunay relative to ip if (1) it can be circumscribed by so called 
empty ellipsoid £ defined by an inequality of the form 

V?[x - c] < R 2 , 

with center c6R d and radius R G M>o, (2) £ has no interior Z d -elements, and (3) vertD = 
£ fl Z d . The collection of all such Delaunay polytopes fit together facet-to-facet to tile M. d , a 
tiling that is uniquely determined by ip. This is the Delaunay tiling of Z d £g>IR relative to 
the form form ip. If a second form t? has Delaunay tiling D#, and if is GL(n, Z)-equivalent 
to D^, then ip and d are forms of the same L-type. 

The description we give below requires that certain degenerate positive forms be admitted 
into the discussion, namely, those forms ip for which ^etu p is a fully rational subspace of 
Z d ® KL The Delaunay polyhedra for such a form are themselves degenerate - they are 
cylinders with axis &zx^ip and Delaunay polytopes as bases. These cylinders fit together 
to form the (degenerate) Delaunay tiling D^. For example, if p G Z d is primitive, i.e. 
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GCD(pi, . . . ,Pd) = 1, then </?[x] = (p • x) 2 is such a form - the kernel £etK is the solution 
set for </?[x] = 0, and given by p -1 , which is fully rational. The Delaunay tiles are infinite 
slabs, each squeezed in between a pair of hyperplanes p • x —k, p • x =k + 1, k G Z. These 
fit together to tile Z d <g> R. 

Let <£> d be the cone of positive definite forms and semidefinite forms with fully rational 
kernels ($ d is sometimes referred to as the rational closure of the cone of positive definite 
forms). For each form ip G $ d there is a Delaunay tiling Del(ip) for 7L d <g> R. 

Definition 5. If D is a Delaunay tiling for Z d , then the following cone of positive definite 
quadratic forms 

$ B = { up e $ d | Del(cp) = V } 

is called an L-type domain. 

For this definition the Delaunay tilings can be the usual ones, where the tiles are Delaunay 
polytopes - or they could be degenerate Delaunay tilings where the tiles are cylinders with 
a common axis K. 

The relatively open faces of an L-type domain, are L-type domains. If D is a trangulation, 
$d has full dimension ( d 1 2 1 ) ; this is the generic case. If D is not a triangulation, dim $d is 
less than , and $d is a boundary cell of a full-dimensional L-type domain. In the case 
where dim$D = 1, the elements p> G $d are called edge forms - Voronoi showed that such 
domains are exactly extreme rays of full-dimensional L-type domains. 

Let 7r$ be the projection operator from CP d onto its quadratic part. If D G 2) and dim D = 
d, then $£> C 7r$(T^,). Since this containment holds for all Delaunay d-tiles D G D, there 
is the following description of $d in terms of inhomogeneous domains: 

$ D = Pi ^(^), 

where the intersection is over all d- dimensional Delaunay polyhedra in D. (It is also true 
that the intersection can be taken over all Delaunay polyhedra in D.) Since the containment 
$d C 71^(7"^) also holds for all Delaunay tilings T) that contain D, there is the following 
description of 7r$(J"^) in terms of homogeneous domains: 

MK) = U ^ 

D3D 

where the disjoint union is over all Delaunay tilings of Z d ® R that contain D. This holds 
not only for full-dimensional cells of 3*^, but for cells of all dimensions. The last equality 
shows that 7:^(7%) is tiled by L-type domains. It also establishes the following 

Proposition 6. If p G "y d is perfect and, therefore, an edge form for an inhomogeneous 
domain, then 7r$p is an edge form for an L-type domain. 

This result can also be established in a more direct way by appealing to the definition of 
L-type domain: if D is a perfect Delaunay polyhedron, then 7r$(CP^,) is a one-dimensional 
L-type domain. By definition, the elements of -K^(7 d D ) are then edge forms. 
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The converse of theorem does not hold - there are edge forms for L-type domains that 
are not inherited from perfect elements in t S' d . Evidence has accumulated recently that 
the growth of numbers of types of edge forms with dimension is very rapid starting in six 
dimensions (see Dutour and Vallentin, 2003), but the growth of perfect inhomogeneous forms 
is much less rapid - there is some hope that a complete classification can be made through 
dimension nine. 

3 Inhomogeneous domains and L- types 
for arbitrary lattices 

In the previous section we only treated the case of "ground lattice" Z d and tilings of Z d <g> R. 
Often, it is more convenient to work with a ground lattice which is a linear or afline section 
of Z d or even a centering (superlattice) of Z d . A point lattice in a real or rational affine 
space is a set L such that for any x Eh the set of vectors {y — x \ y G L} is a free Z-module. 
Z d is a Z-module, i.e. a "vector" lattice, but the "ends of vectors" of Z d form a point 
lattice in Z d <g> R - it is customary to use the symbol Z d in both senses. If A is a section 
of Z d defined by an affine equation, where the Z d is considered as a Z-module, then the set 
A — A = {z — z' | z, z' G A} is a module; if L is an affine section of point lattice Z d , then L is 
a point lattice. The point lattice setup is more intuitive in the study of quadratic functions 
and Delaunay polytopes, while the "vector" (Z-module) setup is more natural for studying 
quadratic forms and, perhaps, Dirichlet-Voronoi cells. A case can be made for either setup 
for Delaunay tilings, for they show both homogeneous and inhomogeneous aspects. We will 
use the term "affine sublattice" in both situations; the distinction will be clear from the 
context or notation. 

For example, one convenient way to describe the root lattice A n is to start from Z n+1 
and consider a sublattice of Z n+1 defined by equation z 1 + ■ — h z n+1 = 1 with respect to the 
Euclidean lattice norm || || (sum of squares) on Z n+1 . Denote by A the lattice of all vectors 
z — z' where both z, z' lie in L = {x G Z n+1 | Yli +1 x i = !}• The resulting pair [A, || | ] is 

/ n+l 

known as the root lattice A n . It can also be regarded as a point lattice I L, ^2 x i k 

A pair [A, tp], where A is a free Z-module of finite rank and tp is a quadratic form on 
A is called a quadratic lattice ( quadratic module is also used). A pair (L, /), where L is 
a point lattice of finite rank and / is a quadratic function on L is called a quadratic point 
lattice. Quadratic lattices [A, tp] and [Ai,</?i] are called isomorphic if there is a Z-module 
isomorphism m : A — > Ai such that </?[x] = </?i[mx] for any x G A. Quadratic point lattices 
(L, /) and (Li, f\) are called isomorphic if there is a point lattice isomorphism m : L — > Li 
such that 7r$/[x] = 7r$/i[mx] for any x G L. Although L does not have a preferred origin, 
the quadratic part of /|l is independent of the choice of such. The notion of Delanay tiling 
makes sense for both types of lattices, the only difference is that in the case of quadratic 
"vector" lattices it has a preferred origin. In the former case we will denote the Delaunay 
tiling of A with respect to tp by Del(A,tp), and in the latter case we denote the Delaunay 
tiling of L with respect to the quadratic part of / by Del(L, 7r$/); when tp or / are defined 
on a superset of a lattice we may omit the explicit respriction in the notation, e.g., we may 
use tp instead of more proper n<s>(p\i^. 
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4 Structure of Perfect Delaunay Polyhedra 

In this section we will describe how to construct new arithmetic types of perfect Delaunay 
polytopes in dimension d from known types of perfect Delaunay polytopes in dimension d — 1 
for d > 2. Since there are no perfect Delaunay d-polytopes for 1 < d < 6, this construction 
can be used starting from d = 7. 

Below tp denotes any positive quadratic form, i.e. tp[z] > for any z. Recall that we 
denote by Tiq>p the quadratic form part of a polynomial p. Our starting point is the following 
structural characterization. 

Theorem 7. (Erdahl, 1992) A polyhedron P G Del(A,tp) is perfect if and only if there is a 
quadratic function p on A with n^p = (p such that Pfl A = V(p) = {v + z | v G vert D, z G T}, 
where D is a perfect polytope in Del (A fl aff D, 7r$(<£>| a g- #)) and r is a submodule of A such 
that A is the direct sum of modules (A D aff D) — (A fl aff P) and T f or, equivalently, point 
lattice A is the direct affine sum of point set A fl aff D and module Y). 

As an example, consider primitive vectors a, b G A = Z d such that ab=l. Let 
</?[z] = (a • z) 2 and let us use a -1 for {x€l d |x-a=0}. Then £ — conv{0, b} is an empty 
perfect ellipsoid in the 1-lattice (aff £q) ^ The set {0, b} = £ fl Z d is the vertex set for a 
perfect Delaunay polytope D in aff(£onZ d ), and V = (a i nZ rf ). We have Z d = Zb©(a ± nZ d ). 
Thus, the polyhedron P is Euclidean direct product of conv{0, b} and a^. The perfect 
quadric circumscribed about P is given, up to a scaling factor, by the equation p(x) = 0, 
where p(x) = (a ■ x)(a ■ x— 1) = </?[x] — a • x. The preimages of negative real values of p 
lie between two hyperplanes a • x =0 and a • x =1, a region whose closure is a degenerate 
perfect Delaunay polyhedron P. 

By the above theorem, for each perfect P G Del(A, ip) the set vertP of lattice points on 
P can be written as vertP©T, where D is a Delaunay polytope and T is a (vector) sublattice 
of A. The direct sum © sign means that (i) each vertex of P can be represented as a vertex 
of D plus a vector of T and (ii) for any fixed choice of the pair (D, T) such representation is 
unique. We call this a direct decomposition of P. If P is a polytope, it is obviously unique. 
What are the possibilities in the case of a degenerate PI 

Lemma 4.1. Suppose vertP = vert-D © V = vertP/ © T are direct decompositions of a 
perfect polyhedron P in Del(A,ip), where D and D' are Delaunay polytopes and Y is a 
(vector) sublattice of A. If&ED is parallel to aff D' , then D' is a translate of D by a vector 
ofT. 

Proof. Since vertP = vertP © T, any v' G vertP/ can be written in a unique way as v + u, 
with v G vertP and u G T. Since aff D' is parallel to aff D, we have aff D' = aff D + u and 
D + u C aff D'. In particular, vert D + u C aff D' n A. So, vert D + uC vert P n aff D' . Since 
vert P = vert D' © T, we have vert P D aff D' = vert P/ + = vert D' and vert D + u C vert P/. 
Since P is a perfect Delaunay polytope in Del(aS D fl A, 7r$((/?| aff £>)), its translate D + u is 
also a perfect Delaunay polytope in Del(aff D' fl A, 7r$(v?| a fr d>))- But this means P + u is 
maximal in Del (aft D' fl A, ir<s>(<p\aKD')) and P' = P + u. □ 

The following theorem gives a universal construction of a perfect Delaunay polytope in 
dimension d + 1 from a perfect Delaunay polytope in dimension d. 
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Theorem 8. Let P be a perfect polytope in Del(A,ip) and let p be its perfect quadratic 
function, i.e. vertP = V(p) and (p = n^p. Suppose D G Del(A, (p) is another Delaunay cell 
of full dimension, which is not a A-translate of P. If e ^ A, then there is a positive definite 
form ip on A © Ze and a perfect polytope P' in Del(A © Ze, ip) such that P' fl aff A = P and 
P'n{affA + e} =D + e 

Proof. We can safely assume G vertP and G vert P. Let then f(y) = ip[y] + /(y) be a 
quadratic function with linear part / such that: (i) /(y) = is the equation of a quadric 
passing through the vertices of P and D + e, (ii) f\ A — p. Since any y G A © Ze can be 
uniquely written as x + ke with x G A and k G Z, we have ^[y] = ^[x] + t/j[ke] +2^(x, ke) = 
<p[x\ + k 2 ifj[e] + 2kip('x, e) and /(y) = /(x) + A;/(e). Z | a is known from p. So, / is determined 
uniquely if and only if ip[e], ip(e, e±), . . . , ip(e, e^), /(e) are fixed. 

We will first show that for any positive (> 0) value of ip[e], function /, that satisfies 
the conditions (i), (ii) is unique, and /(y) < defines a quadric, circumscribed about 
PLi{D + e}, such that no points of A©Ze satisfy f(y) < . Let us fix some ^[e] > 0. Since 
dim D — d, it has d + 1 affinely independent vertices. For each non-zero v G vert D consider 
equation /(v + e) = 0. If Vi, . . . , are independent, then the resulting system of equations 
is independent. (The origin has already been used to set the constant term of / to 0.) We 
know the line joining the circumcenters of D + e and P is perpendicular to A with respect 
to ip. This condition can be written as ip{co — (c P + e), v) = 0, where v G vertL>\0, and 
cd, cp are the circumcenters of D and P. The resulting system of d + 1 linear equations has 
full rank. It is inhomogeneous because </?[v] +V ; [ e ] > 0. So, the solution is unique. Positivity 
of ip follows from positivity of ip, ip[e] > and elementary geometric considerations. For 
a sufficiently small value of ^[e] > the interior of {z G (A © Ze) © H. | /(z) < 0} is free 
of lattice points. This is easy to establish by a calculation, or just by observing that large 
values of ^[e] correspond to ellpsoids that lie inside of an infinite slab of (A©Ze) ©M, which 
is squeezed between hyperplanes corresponding to k — — 1 and k = 2. 

Thus, under the assumptions of the theorem, for a fixed ^[e] > 0, the quadric circum- 
scribed about P U {D + e} is uniquely defined, has positive quadratic part, and is Delaunay 
for sufficiently large -?/>[e]'s. We can now start continuously decreasing ^[e], at a constant 
rate, until the ellipsoid hits a new lattice point, i.e. a point of A © Ze \ (P U {D + e}). 
If it does hit a new lattice element a in a finite time, we are done, since /(a) = will give 
us a linear equation independent of those that we used to determine ^(e, ei), . . . , ^(e, e^) 
and /(e). This new equation will determine ip[e\; the resulting Delaunay ellipsoid /(z) = 
will be perfect. Note that /(z) =0 may contain other new lattice elements besides a. If 
such collision does not happen in finite time, then ^[e] =0 determines a perfect Delaunay 
surface, which is a cylinder with elliptic base. By the preceeding lemma any section of this 
cylinder by a sublattice parallel to A is the vertex set of a lattice translate of P (by a vector 
ne for some n G Z). But the section of A © Ze given by k = 1 is a translate of D, which 
is not a translate of P. It only means that the infinite cylinder scenario is impossible and 
there exists a required /. □ 

Since P is perfect in Del(A,ip), it determines ip (up to a scale factor) and the corre- 
sponding Delaunay tiling of A uniquely. However, there can be different choices of D. For 
example, the 35-tope of Erdahl and Rybnikov lives in a 7-dimensional Delaunay tiling which 
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has many inequivalent (even combinatorially!) types of Delaunay polytopes of full dimen- 
sion. Applying this theorem with different choices of D results, in general but not always, 
in arithmetically distinct perfect polytopes in the next dimension. For example, the De- 
launay tiling determined by P = Gj (this is the Delaunay tiling of [E 7 , || ||]) consists of 
translates of a single copy of GV, and many (more than two - one always has at least 2 
for a simplex) translation classes of Delaunay simplices of index 2 (see Erdahl 1992 for an 
exact description). All these simplices are arithmetically equivalent, in particular, they are 
isometric relative to the form determined by Gj; however, there is more than one orbit of 
such simplicies with respect to the subgroup generated by lattice translations and inversion 
with respect to the origin. But no matter which copy of such a simplex we use for D, the 
resulting perfect Delaunay polytope and quadratic lattice in dimension 8 is the same. On the 
other hand, different choices of translationally-inequivalent -D's in the case of the Delaunay 
tiling determined by the 35-tope sometimes lead to inequivalent perfect Delaunay polytopes 
in dimension 8. 

So, this construction always works, but is not always uniquely determined by P. The only 
known cases where one cannot find a D which is not a translate of P, are those of d = and 
d — 1. A surprising fact is that for d = the construction works even without such a D, since 
[0, 1] is perfect Delaunay in Z d+l for d = 0. We do not know of any other example of a perfect 
Delaunay polytope which can be translated by a lattice vector to every Delaunay polytope 
of full dimension in the Delaunay tiling! Very likely that such an example does not exist. If 
we drop the perfection condition, the only example that we know of is a (^-parallelepiped. 
There should be a way to show that no other examples exist, but, embarassingly, we do not 
have a proof at the moment. 

5 Infinite Sequences of 

Perfect Delaunay polytopes 

To simplify notation we will denote an afline sublattice spanned by all afline integral com- 
binations of a set of lattice points S by aff % S, and we will use the same abbreviation for 
the (vector) lattice spanned by all integral linear combinations of vectors {z — z' | z, z' e S}. 
Denote by j the vector of all ones, i.e. [1, . . . , 1], in Z d . Denote by Vf 2fc the following set in 

vf, 2i M[i*,o-|-^j}xQ, 

where s, k e N = {0, 1, . . .}. All permutations of entries are taken, so that |D^ 2fc | = (j) + 
(s+i) = (^t 1 ) • I n g enera l> whenever we use a shorthand notation for vectors or points such as 
[l s ,0 d_s ], it is understood that all permutations of the components are taken, except for those 
components that are separated from others by semicolumns on both sides, or by a semicolumn 
and a bracket. Also, let us set 

Df, 2fc = { v - v' | v, v G V s % } 

and 
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P* 2k 4 conv { Vf ;2fe U -Vj 2fc U V*. 1>2fc U -Vf +1)2fe } . 
The following is the main constructive theorem of this paper. 

Theorem 9. Let d,s,k G N. If 2 < k and 1 < s < £/ien -Pf 2fc a symmetric perfect 
Delaunay polytope for the quadratic lattice [Af 2 fc, V?s,2fcL where A^ 2fe = aff^ vert P s d 2fc . T/ie 
circumscribing empty ellipsoid is defined as {x G R* 1 1 Vs,2fc[ x ] — -^ 2 }> where 

</?f >2fe [x] = 4k(d-2sk - A;)|x| 2 + (d 2 - (4fc + 2s + l)d + 4A;(2s + fc)) (j • x) 2 . 

The origin is the center of symmetry for P s d fc and does not belong to Af k . Thus, each pair 
of integers (s, fc), for s > 1, k > 2 and s < ^, determines an infinite sequence of symmetric 
perfect Delaunay polytopes, one in each dimension. For s — 1, k — 2 the infinite sequence is 
the one described in the opening commentary, i.e., Gj,d > 7, where the initial term is the 
Gosset polytope G 7 7 = Gj = 3 2 i with 56 vertices. 

The ( d ^ 1 ) diagonal vectors Df 2fc UDf +12fc for Pf 2k have the origin as a common mid-point, 
forming a segment arrangement that generalizes the cross formed by the diagonals of a cross- 
polytope. Moreover, these ( d ^ s 1 ) diagonals are primitive and belong to the same parity class 
for 2fc , namely, they are equivalent modulo 2Af 2k . More generally, primitive vectors u, v 
in some lattice A, with mid-points equivalent modulo A, are necessarily equivalent modulo 
2A and thus belong to the same parity class. And conversely, the mid-points of lattice 
vectors u, v e A belonging to the same parity class are equivalent modulo A. By analogy 
with the case of cross-polytopes, we call any such arrangement of segments or vectors a cross. 
The convex hulls of such crosses often appear as cells in Delaunay tilings - cross polytopes 
are examples, as are the more spectacular symmetric perfect Delaunay polytopes. There 
is a criterion, essentially due to Voronoi (1908), but first formulated and formally proved 
by Baranovskii (1991), that determines whether a cross is Delaunay: Let A be a lattice, <p 
a positive definite form, C the convex hull of a cross of primitive vectors belonging to the 
same parity class. Then C is Delaunay relative to tp if and only if the set of vectors forming 
the cross is the complete set of vectors of minimal length, relative to if, in their parity class. 
This is the criterion we have used to establish the Delaunay property for the symmetric 
perfect Delaunay polytopes Pf 2 k- 

The following result shows that asymmetric perfect Delaunay polytopes can appear as 
sections of symmetric ones. 

Theorem 10. For d > 6 let u = [-1 2 ; l*" 1 ] G Z d+1 . Then 

Gq = conv{ v G vert P^ 1 | v • u = - } 

is an asymmetric perfect Delaunay polytope for (Af 4 , ipf where Mf 4 = afizvertG^. 
The circumscribing empty ellipsoid is defined as {x G affsGg | ^^[x] < R 2 }, where 

<Hx] = 8(d - 5)|x| 2 + (d 2 -9d + 22)0 ' x) 2 . 
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Poly tope 


dimP 


| vertP| 


Symmetry 


pd+i 

1,4 


d+1 




centrally-symmetric 


T d = 


d 


d(d+2) i 
2 1 


asymmetric 



Table 1: Properties of constructed perfect Delaunay polytopes 



[0 d ] x 1 


[-l,0 d - 2 ;l]x(d-l) 


[ld-i._ (d _ 3)]xl 


[0,l*- a ;-(d-4)]x(d-l) 


[l 2 ,0 d - 3 ; 1] x ( rf - 1 K rf - 2 ) 


[l,0 d 2 ;0] x (tZ-1) 



Table 2: Vertices of Gg-topes (same as T d ) generalizing Gosset's G 6 on 27 vertices 



This is the infinite sequence of GVtopes described in the introduction, with Gossett polytope 
Gg = Gq = 22i as the initial term for d = 6. Note that the formula for the form on Mf 4 is 
given with respect to 

The terms in this sequence have similar combinatorial properties. For example, the 
lattice vectors running between vertices all lie on the boundary, in all cases. These lattice 
vectors are either edges of simplicial facets, or diagonals of cross polytope facets-there are 
two types of facets, simplexes and cross polytopes. The Gossett polytope G 6 has 27 5- 
dimensional cross-polytopal facets, but for other members of the series Gg, the number of 
cross-polytopal facets is equal to 2d. Gq can be found as a section of G$, but Gg and Gg do 
not have sections arithmetically equivalent to Gq. 

We summarize the properties of both G-series in Table 1. Table 2 (note that entries 
separated by semicolumns are fixed) gives coordinates of the vertices of Gg-topes (named 
T d -polytopes there), discovered in 2001 by Erdahl and Rybnikov. 

5.1 Delaunay Property Part of Main Theorem 

Let e 1; . . . , stand for the canonical basis of Z d . We will use the following notation: y?i[x] = 
[Yli=i x i) j V 9 2[x] = x— - j ; v 9 2[x] is the squared Euclidean distance from x to the 
line Rj. The following proves the Delaunay property of polytopes P d 2 k asserted by Theorem 

m 

Theorem 11. Let d,s,k E N, and let 2 < k and 1 < s < jr. Let also 

'd> 



^{[(-l^l^x 



A = Z( ei, ...,e d , J 



f ' d-2k, 
Ai = {z G A | z • L 1 = 1 mod 2} 
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Then there is a positive definite quadratic form of the type 



</?[x] = a^i[x] +f3ip 2 [x\, (1) 
with a, (3 G Q>0; such that Pf 2 k ^ s a Delaunay polytope in the lattice (Ai, v?|ai)- 
Below n = d — 2k 

Lemma 5.1. Suppose y?[x] = ay?i[x] + /3</? 2 [ x ] where a, (5 > 0. T/ien a// points zgAi which 
are closest to respect to ip, ie., 

ip[z\ = min{{p(u) : -u G A } 

are, to permutations of components, of the type 

sgn(/)[1^0^'] + ^j, 

with — | < I < | 7 /or some a G Z. Furthermore, each point of A l7 closest to the origin, has 
only one representation of the described type. 

Proof. Suppose z G A 1 is minimal with respect to </?, i.e., </?[z] = miny?. Let z = c^ex + • — h 

Ai 

a d e d + a^, where oi, . . . , o<j, a 6 Z and A = ai + • • • + a^. We have 

d 

E 



y? 2 [z] = y? 2 



a%Q% 



i=i 



a 

^Oid-A- 

i=i 



2 ^ 

42 



Let us prove that the coefficients {<2j} are of two consecutive integer values. Suppose, 
to the contrary, there are and a,- such that — aj > 2. Consider the vector z' = 
cnei + • • • + (oj — l)ej + ••• + (% + l)ej + • • • + a d e d + oA We have y? 2 [z'] - y?2[z] = 
(ai — l) 2 + (a, + l) 2 — af — a 2 = 2(aj — a^) + 2 < — 2 and j • z' = j • z. Since z' G Ai and 
(p[z'] < y?[z], it follows that the vector z is not closest to z which is a contradiction. 

Now, let b be the smallest of the values of the coefficients {a^}. Subtract oj from the first 
part and add an equal value of bk± to the second part of the existing representation of z 
as an integral linear combination of e^'s and ^. After a permutation of the components z is 
equal to sgn(/)[l |z| , d_|i| ] + (a + ofc) 1 where < / < d. If I > \ , again subtract [l d ] from the 
first summand and add j to the second summand to get the required representation. 

Note that although dim A = d, we represent vectors of Ai as integral linear combinations 
of d + 1 vectors; thus, we have to understand if the representation described in the statement 
of the lemma is unique for each of the points of A x that are closest to the origin. To prove 
that each minimal vector of A , with respect to ip, has only one encoding of the described 
type, note that the components of the vector sgn(/)[l''l; O^ - '''] are of at most two values. 
In sgn(/)[l''l; d ~'''], either l's fill the positions of the largest value, or — l's fill the positions 
of the smallest value. Since — | < I < |, the choice is unique. □ 

Lemma 5.2. All points of the affine lattice A±, which are closest to with respect to a form 
ip = atpi + (3^2, where a,(3 > 0, belong to the set {z G Ai | |z • j| < ^}. If z is closest to 
with respect to y? and \z ■ j| = ^, then z G 
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Proof. By multiplying both a and (3 by the same positive number, we can assume that the 
minimal value of cp on Ai is 1. Since - G A 1; (p[-) = a(-) 2 > 1. Let z G IR d be a point with 

tp[z\ = 1. Represent z = + u where 7 G R, u • j = 0. We have tp[z] = a^—^ — V /3|u| 2 = 1, 
therefore a-^- < 1. Since ck(^) 2 > 1, we have proven that I7I < 1, so |z ■ j| = |7|- < -. 

If the latter inequality holds strictly, then | | = 1 and ip[z] = + (3\u\ 2 = 1. Since 
a (^) 2 — 1) we necessarily have (3 = so z = ±-. □ 

Proof of Theorem 6: Let <p a ^ be the type ([I]), i.e. W/3[ x ] = ot^ifx] + /^^[x]. Let us 
first describe a "reduced" set of points M, such that any point of Ai where the minimum of 
<Pa,p over Ai is attained, belongs to this set, possibly after a permutation of the components. 
Then we "handpick" the form (p = (p a ^ so that the vertices of Pf 2 k are the minimal points 
and other points of Ai are not. Below is an implementation of this plan. 

Consider the set 



M 



\zeA 1 z = sgn(/)[l"l; 0^1] + J-, < z • j < -, -~ < I < ±) U {!}. (3) 
[ n n 2 2 J n 

By the two previous lemmas, every z G Ai, which is closest to the origin of A with respect 
to ip, perhaps after a permutation the components, belongs to this set. 
For z = sgn(/)[l^; d ~\ l \] + a-, we have 

ad 

z j = H , 

n 

{ d\ 2 

<pM = ( l + —) 



n J (4) 

y 2 [z] = |Z| 



Z 2 
d 1 

L 1 ■ z = a + I mod 2. 



From the calculations above, we get that M 
M 



{sgn(Z)[l |,| ;0 <H,| l +a'- l + a = l mod 2, < In + ad < d, -- < I < - \ U {^-}. 
[ n 2 2 J n 

(5) 

Thus, we have shown that 

(A) For any point of A 1; minimal with respect to (p a ,p 7 there is a permutation of the 
components that turns it into a member of M given by 

We define a mapping from M. d to M 2 by m(x) = (y?i[x], ^[x]). We will call the image of 
M the (p- diagram. Lines parallel to the <£>i-axis in M. 2 will be called horizontal. Next, we will 
show that 

(B) //zi,z 2 G M, then m(zi) and vniz^) on the ip-diagram belong to the same horizontal 
line if and only if Z\ G ±z 2 . 

To prove this, take z = sgn(Z) [l 1 ' 1 ; d_|/| ] + a± G M. If Z ^ 0, then the condition < 
In + ad < d uniquely determines a for a given value of I, and if Z = 0, then a — 1. Therefore 
Z uniquely defines a. 
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The function I — > |Z| — ^ is even and increasing on [0, |], so two different points of the 
93-diagram may belong to the same horizontal line if and only if their preimages are zi = 
[l |,| ,0 d-|,| ]+aii and z 2 = -[l l ' l ,0 d_|i| ] + a 2 i for some l,a u a 2 . If I ^ 0, then < lk + a x d < d, 

< —lk + a 2 d < d. Adding these inequalities, we get < (ai + a 2 )d < 2d, so a\ + a 2 G {0, 1}. 
If ai + a 2 = 0, then Zi = — z 2 so m(zi) = m(z 2 ). If ai + a 2 = 1, then the numbers I + a\ and 
—I + a 2 have different parity which contradicts conditions I + a x = —I + a 2 = \ mod 2. If 

1 = 0, then ai = a 2 = 1. This proves claim (B). 

(C) The proof of Delaunay property is based on that each line ax\ + f3x 2 = 1, where 
a, (3 > 0, which contains an edge of the convex hull of the 92-diagram, gives rise to a quadratic 
form </?[x] = cc</?i[x] + /3</? 2 [x] such that the ellipsoid tp[x\ < 1 passes through points of the 
affine lattice Ai and does not have any points of Ai in its interior. 

We are going to show that the following set is the vertex set for a centrally-symmetric 
Delaunay polytope in Ai with center at the origin of A. 

Vf, 2fc U V d s+1 , 2k U -Vj 2fc U -V d s+ i,2k (6) 

A point from Vf 2fc will be denoted by v Sj2 fc. 

(D) First consider the case of d — 7, 2k — 4, s — 1. The (/9-diagram for these values of 
parameters is shown in Figure 1 (left). We see that the line |xi + \x 2 = 1 passes through 
points m(vi 5 4) = m([l,0 6 ]) and m(v 2) 4) = m([— 1 2 ,0 5 ] + |), and all other points of the tp- 
diagram are contained in the open half-plane |xi + |x 2 > 1. This means that polytope P/ 4 
is a Delaunay polytope with respect to quadratic form fy?i[x] + |<^ 2 [x] = |x ■ x + |(j ■ x) 2 . 

(E) Now suppose that d > 8, 2 < k, and 1 < s < 4=£=±. Suppose z = sgn(Z) [l' 1 ', O^l'l] + 
G M and < (p 2 [z] < ^ — (or, equivalently, |/| < ^). We want to prove the following 

two statements (El) and (E2). 

(El) If\l\ = 7^ (clearly, in this case ^ is an integer), then 

± .= ± ([l*,0^*]-(|-l)i). (7) 

A direct check shows these points belong to set M. Then we use claim (B) that there are at 
most two points that can be written as 

[lKiy-KI] +a i e M 
n 

for any fixed value of |/|, and there are excatly two if and only if Ik + ad = 0, in which case 
the points are ±z. 

(E2) If0<l<4:, then 

z= [il<l,0 d -KI]-(Z-i)i (8) 
Recall that v Sj2fc stands for [l s ,0 d ~ s ] — (s — 1)£. We have 

/ d\ 2 s 2 

Vl[v s ,2fc] = ( S + (1 - S)- J , ^2[v s , 2 fc] = S - — . (9) 
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Figure 1: ^-diagrams for d = 7 (left) and d = 8 (right); 2k = 4 



All points m(v S 2fc) therefore belong to the parabola 

'"(('^-"S''-?) (io) 

which touches the vertical axis in the point (0, ^ — ^p) when £ = The parabola is shown 
in dash on Figure 1. The portion of the parabola for < t < is the graph of a convex 
function. 

Therefore, for each 1 < s < ■M- we can find a line with equation axi+/3x 2 = 1 which passes 
through points m(v Sj2fc ), m ( v s+i,2fc), supports the convex hull of the (^-diagram and does not 
contain points m(z) for z e M \ {±v St2 k U ±v s+ i j2 fc}. Quadratic form <p[x] = a^i[x] + /3<£>2[x] 
defines an empty ellipsoid centered at 0, which contains the vertices of Pf 2 k on ^ s boundary, 
and does not contain any other points of Aj. 

We have proven that -P,f 2fc * s a Delaunay polytope in affine lattice Ai with respect to 
quadratic form if = ^ M for d > 7. Explicit formula (jHJ) for 2 fc is established by a direct 
calculation. 

An example of the yj-diagram for d = 19, 2k = 6 is shown in Figure 2 (right-hand image). 
Note that not all points of m(M rec ) belong to the parabola (TTUT) . 
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Figure 2: ^-diagrams for d = 9 (left) and d — 19 (right); 2k = 4 



5.2 Perfection Property Part of Main Theorem 

Theorem 12. Let d, k, n = d — k, s, a<pi +/?v 9 2 as in Theorem^ Then, there is at most 
one pair consisting of a quadratic form <p = atpi + fip>i and c G M d such that the quadric 
(p\x. — c] = 1 circumscribes the polytope Pf 2 k ■ 

Proof. Since dim P/ 2fc = d and -P,f 2fc nas as the center of symmetry, is the center of sym- 
metry of the circumscribing quadric. Thus, c = 0. Therefore, if <p and if) are two quadratic 
forms such that the corresponding quadrics circumscribe -P s d 2fc , then (ip — V ; ) I v d 2fc uv d +1 2k = 0- 
Let us consider an arbitrary quadratic form £ such that £lv d 2fc uv d +1 2k = an d prove that 
£ = 0: this will prove the theorem. 

We will use the following symmetrization technique. Let G be a subgroup of the group Sd 
of permutations on the d coordinates. If r £ G, then we write t : i *—>■ ir. The symmetrization 
of the form £ by G is defined as 



Since the polytope P s d 2fc C M. d is invariant under permutations of the coordinates of M. d , all 
components in the above sum are on \iertPf 2k . Therefore £lv d 2fc uv d +1 2k = 0- 

First we prove that if a form £ satisfies the condition £|v d 2fc uv d 2k = 0> then the sum of 
the diagonal elements of /, and the sum of off-diagonal elements are 0. We symmetrize / by 
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the group G = Sd and get a form Sym G £ with diagonal coefficients proportional to the sum 
of diagonal coefficients of £, and non-diagonal coefficients proportional to the sum of non- 
diagonal coefficients of /. Therefore Sym G £[x] can be written as a</?i[x] + /3</?2[x]. Suppose 
that a and (3 are not both equal to 0, i.e. a 2 + f3 2 > 0. Since Sym G £| V d uV d = 0, 

the line axi + j3x 2 = passes through the points m(V^ 2lt ) and m{Vg +12k ), where m(x) = 
(</?i[x], (/5 2 [x]). This is impossible because the line that passes through these points is uniquely 
defined and does not contain 0. This contradiction proves our claim. 

Next we prove that all diagonal elements of the form £ are equal to 0. Suppose that one 
of them is nonzero. Without a limitation of generality we may assume that t — fu ^ 0. Let 
St(l) be the subgroup of permutations which leave 1 fixed. The matrix of Sym^) £ is 



t 




(3 




■ P 


p 


a 


5 




. 5 


p 


5 


a 


5 


. 5 




5 




5 


a 5 




5 






5 a 



(12) 



Consider the following vectors: 

Vl = [1-; d -°] - (s - l)i G V d s , 2k v 2 = O^ 1 ] - 8 ± E V*. 1>2fc . 

The conditions Sym St ^ £(vi) = 0, Sym St ^ £(v 2 ) = yield 

t + 2(s- 1)13 + (s - l)a + (s- l)(s - 2)5- 

2(s -l)(t + (s + d-2)P + (s- l)a + (s - l)(d - 2)5) 

n 

(s - l) 2 (t + 2{d - 1)13 + (d- l)a + (d-l)(d- 2)5) 

n 2 



, n a , 2s(t + (s + d-l)P + sa + s(d-2)5) 
t + 2s(3 + sa + s(s - 1)5 — — + 

n 



s 2 (t + 2(d - 1)0 +{d- l)q + (d-l)(d- 2)5) = Q 



We also have the conditions that the sums of the diagonal and the off-diagonal elements 
are equal to 0: 

t + (d-l)a = 0, (d -l)(d- 2)5 + 2(d - l)f3 = 0. (15) 
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The determinant of the system of these 4 equations in variables a, (3, 5 and t is equal to 



-(d - - d + 1)0 - d)(d - 2 - n) 



n 



(16) 



and it is not equal to for the specified values of d, n and s. Hence t = 0. 

We have shown that the diagonal coefficients of £ are all zero. Next we prove that all 
off-diagonal elements of £ are also equal to 0. Without loss of generality we can assume that 
£12 7^ 0. Let Synw-n £ be the symmetrization of £ by the group of permutations which map 
the set {1, 2} onto itself. The matrix of Synw-n £ is then 



a (3 
a (3 
(3 p 
5 



P 
(3 (3 



5 



P 
5 





5 



8 




where a^0. We consider the following vectors: 



n 



S ~ ^ ^"s+l,2fc 

n 



The equations Sym^^ £[ui] = 0, Sym St ^ £[u 2 ] = yield 



St(l) 



(17) 



4(s - l)/3 + 2a + - 1)0 - 2 ) 5 ~ 2 -( 2 « + 2(s - 1 + d - 2)/3 + (s - - 3)5)+ 



n 



- l) 2 (2a + 4(d - 2)/3 + (d - 2)(d - 3)5) 



fa + 1)fl(5 _ 2( a + l) a (2/?+(rf-3)fl + s 2 (2q + 4(^-2)1+^-2)^ -3)5) = Q ^ 



n 



We also know that the sum of off-diagonal elements of the matrix of Synw^ £ is equal to 0: 

2a + 4(d-2)/3+(d-2)(d- 3)5 = (20) 
so, with some simplifications, the previous two equations can be rewritten as 



4(s - l)f3 + 2a + - 1)0 - 2)5 - 2-(2a + 2(s + d- 3)/3 + (s-l)(d- 3)5) = 0, 

n 

n 

The systems of equations (1201 and (j2"T|) in variables a, /3 and 5 has determinant 8(d — 
2 — n)(—d + s + 1) which is not equal to for the specified parameters d, n and s which 
proves that the system has only zero solution. In particular, a = 0. We have proven that all 
off-diagonal elements of the matrix of form £ are equal to 0. □ 
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